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We apply the methods of Chiral Perturbation Theory to the analysis of the first moments of the 
Generalized Parton Distributions in a Nucleon, usually known as generalized form factors. These 
quantities are currently also under investigation in Lattice QCD analyses of baryon structure, 
providing simulation results at large quark masses to be extrapolated to the "real world" via Chi- 
ral Effective Field Theory. We have performed a leading-one-loop calculation in the covariant 
framework of Baryon Chiral Perturbation Theory (BChPT), predicting both the momentum and 
the quark-mass dependence for all the vector and axial (generalized) form factors. In particular 
we discuss the results for the limit of vanishing four-momentum transfer where the GPD-moments 
reduce to the well known moments of Parton Distribution Functions (PDFs). We fit our results to 
available lattice QCD data, extrapolating down to the physical point. We conclude by presenting 
outstanding results from a combined fit to different GPDs-moments. 
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1. Introduction 

In this talk we discuss the findings of references [P and where the generalized form factors 
of the nucleon were analyzed in the framework of covariant Baryon Chiral Perturbation Theory 
(BChPT). In standard SU(2) BChPT the results for form factors typically depend on two variables, 
the momentum transfer squared and the quark-mass and on a number of low energy constants 
(LECs). In this work we study the quark-mass dependence of the isoscalar- and isovector-vector as 
well as the isovector-axial generalized form factor A2,o(?) at zero momentum transfer and predict 
the physical value of those quantities by determining previously unknown LECs by a fit of the 
BChPT results to lattice QCD data. 

Working in twist-2 approximation, the parity-even part of the structure of the nucleon is en- 
coded via two Generalized Parton Distribution functions (GPDs) //^(jc, B, , t) and E'^ {x,^,t) 

Moments of GPDs can be interpreted much easier and are connected to well-established 
hadron structure observables. E.g. the zero-th order (Mellin-) moments in the variable x corre- 
spond to the contribution of quark q to the well known Dirac and Pauli form factors Fi (t), F2{t) of 
the nucleon: ^ 

j dxx^m{x,£,,t)=F^{t) , 1 dx:PE'i{x,£,,t)=F^{t). (1.1) 

Our aim is the application of the methods of ChPT to the analysis of the. first moments in x of these 
nucleon GPDs 

fdxxm{x,^,t) =A|o(0 + (-2^)'cf,o(0 , f^dxxEHx,^,t)=Bl,{t)-{-2^fc',,,{t) 

(1.2) 

where one encounters three generalized form factors A2Q{t), 82^(1), C2o(0 of the nucleon for 
each quark flavor q. For the case of 2 light flavours the generalized isoscalar (u+d) and isovector 
(u-d) form factors have been already studied in a series of papers at leading-one-loop order in the 
non-relativistic framework of Heavy Baryon ChPT (HBChPT) We will provide the first anal- 
ysis of these generalized form factors utilizing the methods of covariant BChPT for 2 light flavors 

Our BChPT formalism ^ makes use of a variant of Infrared Regularization [0] for the loop 
diagrams and is constructed in such a way that we exactly reproduce the corresponding HBChPT 
result of the same chiral order in the hmit of small pion masses. For the complete analytical expres- 
sions of the discussed results and a detailed description of the formalism used for the calculation 
we refer to [|l|] and [Q]. 



2. The Generalized Form Factors of the Nucleon in ChPT 



In ChPT one can direcly access the isoscalar (s) and isovector (v) contribution to the general- 
ized form factors of the nucleon by evaluating the following matrix elements [^] : 



i{p'\lY{nD v}q\p)u+d = u{p) 
i{p'\qY{ii'^ v}q\p)u-d = uip') 
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Figure 1: Loop diagrams contributing to the first moments of the GPDs of a nucleon at leading-one-loop 
order in BChPT. The solid and dashed lines represent nucleon and pion propagators respectively. The solid 
dot denotes a coupling to an external tensor field. 



The brackets {. . .} denote the completely symmetrized and traceless combination of all indices 
in an operator, u (u) is a Dirac spinor of the incoming (outgoing) nucleon of mass M^, for which 
the quark matrix-element is evaluated. 

From a powercounting analysis we find that the Feynman diagrams contributing to the first 
moments of GPDs of a nucleon at leading-one-loop order [^{p^)] in covariant ChPT are the ones 
depicted in Fig.|l]. 

3. Analysis of the results for q ^ {t = 0) 

In this section we present an analysis of the generalized isovector- and isoscalar-vector form 
factors ()(?) in the forward limit f — > 0. The details of the ChPT calculation as well as a complete 
analysis of the form factors BY^it), CYoit) and their connection to the spin physics sector can be 
found in ref[[T]]. 

In the forward limit f — > the generalized form factors q {t = 0) can be understood as mo- 
ments of the ordinary Parton Distribution Functions (PDFs) q{x), q{x) 

{x)u±d = A'ilit = 0) = / dxx{q{x)+q{x))^^a . (3.1) 
JO 

Experimental results are available for (x) in proton- and "neutron-" targets, from which one can 
estimate the isoscalar and isovector quark contributions at the physical point [Q] at a regularization 
scale pL. We choose /i = 2 GeV for our comparisons with phenomenology. 
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Figure 2: Left panel. Fitlofthe j.ggjjjjQf gq (3 2)totheLHPClatticedataofref.|[To|. Right panel. 

Fit II of the ^(/:>2) BChPT result of Eq.(3.2) to the LHPC lattice data of ref.[|l§ and to the physical point 
(solid line). The dashed curve shown corresponds to ff{p^) result in the HBChPT truncation. 





^2.0 




4(lGeV) 


Fit I (4 points - 2 parameter) 
Fit II (6+1 points - 3 parameter) 


0.157 ±0.006 
0.141 ±0.0057 


0.210 (fixed) 
0.144 ±0.034 


-0.283 ±0.011 
-0.213 ± 0.03 



Table 1: Values of the couplings resulting from the two fits to the LHPC lattice data for (x)„_^ The 
errors shown are only statistical and do neither include uncertainties from possible higher order corrections 
in ChEFT nor from systematic uncertainties connected with the lattice simulation. 



For the PDF-moment o(f = 0) we obtain to 



>2) in BChPT 



tJ- ( 2 2 / 2 ^ 

9 X fnt. 9 9 mi / mt, 



— '^2,0 ■ 



-8^±^ arccos — ^ 



4fl + 31og4^ 



14-8^ 
Ml 

mi ml 2mj,{AMl - ml)^- ( m^^ 

J log — 5- H -. arccos — — 

M4 ^Ml Ml \2Mq, 



+4m^i^ 
+ ^ Ml 



(3.2) 



Most LECs in this expression are well known from analyses of chiral extrapolation functions 
[^]. However, the sizes of c ^'^2 ^^'^ '^s ' (■^) ^^^^y Poorly known at this point. The coupling 
Afljo is related to the spin-dependent analogue of the mean momentum fraction, namely (Ax)„_rf 
(see Section 4 and [^). In a first fit (Fit I) to lattice data we constrain Aajo from the phenomeno- 
logical value of {Ax)'^'"'"f r::0.21 and perform a 2-parameter fit with the couplings (2205 Cg '(IGeV) 
at the regularization scale A = 1 GeV. We fit to the LHPC lattice data for this quantity as given 
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in ref.[10], including lattice data up to effective pion masses of niji ^ 600 MeV. The resulting 
values for the fit parameters together with their statistical errors are given in table |I] and the re- 
sulting chiral extrapolation function is shown as the solid line in the left hand side of Fig.§. The 
extrapolation curve tends towards smaller values for small quark-masses, but does not quite reach 
the phenomenological value at the physical point, which is not included in the fit. Since BChPT 
is based on a systematic perturbative expansion, it does not only provide us with the result at a 
certain order, but also allows for an estimate of possible higher order effects. From dimensional 
analysis we know that the leading chiral contribution to (x)„_d beyond our calculation takes the 

3 

form ^(p ) ~ — h Constraining between values — 1, . . . , -t-1 (the natural scale of all 

couplings in the observables considered here is below 1) and repeating the fit with this uncertainty 
term included leads to the grey band indicated in FigJ^ As one can see the phenomenological 
value for {x)u-d lies well within that band of possible next-order corrections, giving us no indica- 
tion that something may be inconsistent with the large values for typically found in lattice 
QCD simulations for large quark-masses. 

As stated in the Introduction, the covariant BChPT scheme used in this analysis is able to repro- 
duce exactly the corresponding non-relativistic HBChPT result at the same order by the appropriate 
truncation in 1 / {Idlf' F^Mq) . In order to also compare the ff{p^) HBChPT result of refs.ljn]] with 
the ^{p^) covariant BChPT result of Eq.(|3^) we perform a second fit (Fit II): We fit the covariant 



expression for {x)u-d of Eq.(3.2) again to the LHPC lattice data and we constrain the coupling 
ii such a way, that the resulting chiral extrapolation curve reproduces the phenomenological 
value of {x)f_!^- = 0.160 ±0.006 % exacdy for physical quark masses. The parameter values for 
this Fit II are again given in table [l], whereas the resulting chiral extrapolation is shown as the solid 
line in the right hand side of Fig.^. We would like to emphasize that the curve looks very rea- 
sonable, connecting the physical point with the lattice data of the LHPC collaboration in a smooth 
fashion. For the comparison with HBChPT we now utilize the very same values for ^2 o ^'^^ ^8 
Fit II. The resulting curve based on the &{p^) HBChPT truncation is shown as the dashed curve 
in Fig.^ One observes that this leading-one-loop HBChPT expression agrees with the covariant 
result between the chiral limit and the physical point, but is not able to extrapolate on towards the 
lattice data. 

We therefore conclude that the smooth extrapolation behaviour of the covariant &{p^) BChPT 



expression for {x)u-d of Eq.(3.2) between the chiral limit and the region of present lattice QCD 



data is due to an infinite tower of j terms. 

The same analysis can also be performed for the isoscalar generalized form factor A\Q{t) 
which in the forward limit reduces to q(0) = {x)u+d- Fig.|| shows a 2-parameter fit of the ^{p^) 
covariant BChPT results for this observable [|l|] to LHPC and QCDSF data of reference [ [To| ] and 



[12]. We want to stress that the experimental value of the quantity (x)„+^ is not included in the fit. 
As the plot shows, the obtained chiral extrapolation curve is very satisfying, consistently linking the 
lattice data at large quark-masses with the phenomenological value. In contrast, the correspondent 
result in the heavy baryon limit represented by the dashed line in Fig.^ not even allows for an 
interpolation between the presently available lattice data and results from experiments. 
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Figure 3: Two-parameters Fit of the ff{p^) BChPT result of ref.[0| to the LHPC lattice data of ref.Q and 
to the QCDSF data of ref.[|l2ll. We obtain a^^ = 0.513 ± 0.006 and eg = -0.064 ± 0.005 as values for the 
free parameters. The band shown indicate estimate of higher order possible corrections. The dashed line 
correspond to the respective HBChPT results at this order 

4. Combined Fit 

We have extended the analysis to the first moments of the axial GPDs H'i{x, ,t) and E'^ (x, ^ , t) 

dxxH^x,^,t) =Af o(0 , f dxxE'i{x,^,t)=B%{t). (4.1) 

Again, in the limit of vanishing four-momentum transfer the isovector form factor A\^^{t 0) is 
directly connected to the spin dependent analogue of the mean momentum fraction {Ax)u-d 

^2,o(0) = {Ax)u-d = Jodxx{qi{x)-q^{x))l_^ (4.2) 

where Aaj q corresponds to the chiral limit value of (Ax)„_j. 

Looking at the ^{p^) BChPT expression for {Ax)u-d [Hi one can easily observe that each 
isovector moment ((x)„_j and (Ax)„_d) depends on 3 unknown parameters: 2 couplings (a2 0' 

q) and one counterterm. As the same couplings contribute in both moments, it is hoped that 
a simultaneous fit of our BChPT results to the lattice data of ref . [[I0|] can considerably reduce the 
statistical errors. As one can see from Fig.^, the results of this procedure are pretty outstanding, 
given that the values at the physical pion mass were not included in the fit! The chiral curvature 
in both observables naturally bends down to the phenomenological value for lighter quark masses, 
leading to a very satisfactory extrapolation curve. 

We would like to stress that this efficient cross-talk between the ChPT results for {x)u-d and 
{Ax)u-d occurs only in the covariant framework, while in the non-relativistic approach both ob- 
servables are completely independent at this order. 

We conclude that combined fits of several observables characterized by a common subset of 
ChEFT couplings are the winning strategy towards the most reliable chiral extrapolations of lattice 
QCD results. 
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mJGeV] mJGeV] 

Figure 4: Combined FIT of the &{p^) results of ref. [§] to the lattice data of ref. @. Note that the 
phenomenological values at physical pion mass were not included in the fit. The bands shown indicate 
estimate of higher order possible corrections. 
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